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O ' Abstract 

^ \ Hadron production in leading order pQCD is reviewed. The shape of the single 

^ \ inclusive particle spectra is well described for > 2 — 3 GeV at center of 

Q ■ mass energies from 20 GeV to 2 TeV. The phenomenological K-factor is found 

to decrease systematically with y/s. For ultra-relativistic heavy ion reactions 
the calculation is augmented with the effects of initial multiple parton scatter- 
I ■ ing and final state radiative energy loss. Baseline CERN-LHC predictions for 

^ . hadron production m p + p and suppression in central Pb + Ph reactions at 

^ \ y/s = 5.5 TeV are given in comparison to the corresponding results at BNL- 

^ ■ RHIC and CERN-SPS energies. 

(N : 
^ , 

(N 1. INTRODUCTION 

o ■ 

' One of the main goals of the upcoming p + p program at y/s = 14 TeV at the Large Hadron Collider 

p • (LHC) at CERN is the unambiguous discovery of the Higgs boson, predicted by the standard model of 

' ! particle interactions, as well as the search for physics that reaches beyond our current understanding of 

Qj I the constituents of matter and the force mediators. Equally important, however, is the continuing effort to 

• investigate the strong sector of the SM and probe experimentally some of the fundamental predictions of 

J> , QCD: the deconfinement phase transition and chiral symmetry restoration. To order (a^*^)^ the strong 
coupling constant Ug = gU 47r reads 
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where the lowest order a^^ {Q"^ ^nA = A-k/ (ll — \nf \ In -r^ — . QCD is thus "asymptotically free", 

i.e. for nj < 6 the running coupling, Eq. ([l]), approaches zero in the limit of large momentum transfer 
Q, - a general feature of non-Abelian gauge theories with sufficiently small number of fermions. In hot 
and dense matter the typical momentum scale is on the order of the temperature Q ~ T, assuming local 
thermal equilibrium. Large initial temperatures and energy densities can be experimentally achieved 
in ultra-relativistic heavy ion reactions. The Pb + Pb program at \/sj^j^ = 5.5 TeV at the LHC is 
targeted at the search for the quark-gluon plasma (QGP) ^ and the study of its properties. The Bjorken- 
estimated Tj ~ 1 GeV for these conditions exceeds by a large margin the current lattice QCD results 
for the critical temperature Tc ~ 170 MeV An important advantage of the LHC is that it will ultimately 
make small-x physics studies in heavy ion collisions feasible (not applicable at RHIC except possibly 
for pt < 0.5 GeV). At very small values of x the rapid growth of the gluon distribution in nucleons and 

* Contribution to the CERN Yellow Report on Hard Probes in Heavy Ion Collisions at the LHC. This manuscript also 
provides details on the background calculations for Ref. |Q] and discussion on results that have been omitted there due to space 
limitations. 



nuclei is tamed by absorption terms that lead to a modification [g] of the DGLAP evolution equations 
and correct the small-x unitarity problem: 

dxG{x,Q'^) Cachs dx' x „„ f x 
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where is the gluon splitting function and G^'^\x, Q"^) is proportional to the gluon density overlap. 
An opportunity to test these predictions at the LHC is provided by the p + Pb program at ^/s = 8.8 TeV 
since it stands the best chance of identifying initial state nuclear effects and separating them [Q] from 
the final state multi-parton interactions. 

The search of novel physical effects in ultra-relativistic heavy ion reactions at the LHC can only 
rely on a detailed comparison between the experimental data and the projected current, or "conventional", 
knowledge. This calls for detailed baseline calculations of jet and hadron production at those center of 
mass energies as well as an estimate of the known nuclear effects. The purpose of this manuscript 
is to present a lowest order (LO) analysis of inclusive hadron production up to the Tevatron energies 
and discuss hadron differential cross sections and composition at the LHC. This choice is dictated by the 
requirement of self-consistent incorporation of nuclear effects that are at present computed/parameterized 
to LO. It also complements next-to-leading order (NLO) calculations of jet and hadron production [^, ^. 
Evaluation of the nuclear modification factors at the LHC in comparison to RHIC and discussion of the 
hadron composition is also presented. Results on Cronin and shadowing effects are given at forward (in 
the direction of the proton/deuteron beam) y = +3 rapidity. 



2. HADRON PRODUCTION IN FACTORIZED PQCD 

The standard factorized pQCD hadron production formalism expresses the differential hadron cross 
section inA^ + A^^ /i + Xasa convolution of the measured parton distribution functions (PDFs) 
fa/N{xa, Qa) for the interacting partons (a = a, h), with the fragmentation function (FFs) Df^^ci^, Ql) 
for the leading scattered parton c into a hadron of flavor h and the parton-parton differential cross sections 
for the elementary sub-process da^"'^^'^'^^ / di: 

da^N I } } 

^h—l^ = KnLO^ dZc / dXadXb fa/p{Xa,Ql)fb/p{Xb,Qb) 

abed " 



mill min 



s da^^^^^^^ 

X Dh/c{zc, Ql) — + n + . (3) 

A list of the lowest order partonic cross sections can be found in [^. In Eq. @ Xa,Xb are the initial 
momentum fractions carried by the interacting partons and Zc = Ph/Pc is the momentum fraction of the 
observed hadron. Kj\[ lq is a phenomenological factor that is meant to account for next-to-leading order 
(NLO) corrections. It is ^/s and scale dependent and takes typical values ~ 1 — 4. One usually finds 
that Eq. (^ over-predicts the curvature of the inclusive hadron spectra \dpj,da^\ at transverse momenta 
< 4 GeV. This can be partly corrected by the introduction of a small intrinsic (or primordial) kx- 
smearing of partons, transversely to the collision axis, and generalized parton distributions fa{x, kx, Q"^) 
motivated by the pQCD initial state radiation. For the corresponding modification of the kinematics in 
(^ in addition to the / d^kj, J d'^kj, {■ ■ ■) integrations see [Q]. The generalized parton distributions are 
often approximated as 

/CT, Q^) ~ Q^)g{kT), gikr) = ^ Ji^^2^ > (4) 
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where the width {k^) of the Gaussian enters as a phenomenological parameter. 
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Fig. 1: Extracted KffLo from comparison of LO pQCD calculation to data [ p^ plj , |2| |l3| |T^, |l5|| at and about 
mid-rapidity in the range 2 < pr < 25 GeV. A systematic decrease of Kn lo with \/s is observed and illustrated 
in the bottom right panel. The projected 50% uncertainty at y/s — 5.5 — 8.8 TeV is also shown. 
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Perturbative QCD fits to data [|8|, |[ |T|, [IT], use different coupled choices for Knlo 

and (fcj.) and tlie extracted values are thus not directly comparable. However, similar agreement between 
data and theory at the level of spectral shapes and the ^/s dependence of the corrective factors discussed 
above is found. In [ p^ the factorization and fragmentation scales were set to Qpdf = Pt/2 and 
Qff = Pt/'^Zc and no Knlo factors were employed. The extracted {k'^) decreases from 2.7 GeV^ at 
^/s ~ 50 GeV to 0.75 GeV^ at ~ 2 TeV. Alternatively, in [17] no primordial /c^-smearing was used 
and the scales in the calculation were fixed to be Qpdf = Qff = Pt- The deduced Knlo decreases 
from ~ 6 at - 50 GeV to ~ 1.5 at ^ ~ 2 TeV. 

In the fits shown in Fig. I we have used the GRV98 LO PDFs ^ and the BKK LO FFs ||T^. 
Proton+antiproton fragmentation has been parameterized as in [|^], inspired from PYTHIA [21] results. 
A fixed {k^)pp = 1.8 GeV^ has been employed, leading to a Knlo parameter that naturally exhibits 
a smaller variation with ^/s. A ±25% error band about the Knlo value, fixed by the requirement to 
match the moderate- and high-pr behavior of the data, is also shown. The fragmentation and factorization 
scales were fixed as in p7|]. In the lower right panel the systematic decrease of the next-to-leading order 
K-factor is presented. Two fits to Knlo have been used: linear Knlo = 2.7924 — 0.09991ns and 
quadratic Knlo = 3.8444 — 0.3234 In s + 0.0107 In^ s in In s. For center of mass energies up to 1 TeV 
the two parameterization differ by less than 15% but this difference is seen to grow to 30%-50% at 
./i = 5 - 10 TeV. 




Fig. 2: The predicted LO differential cross section da^^ /dycPpT for inclusive neutral pion and charged hadron 
production at midrapidity y = in p + p (p + p) reactions is shown for y/s — 17, 200, and 5500 GeV. The ratio of 
neutral pions to inclusive charged hadrons versus px is given in the right panel. 



In Fig. ^ the predicted transverse momentum distribution of neutral pions and inclusive charged 
hadrons is shown, corresponding to the quadratic in In s fit to Knlo for energies typical of SPS, RHIC, 
and the LHC. The significant hardening of the spectra with ^/s has two important consequences for p+A 
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and A + A collisions: a notably reduced sensitivity to initial state kinematic effects (smaller Cronin) and 
larger variation of the manifested final-state multi-parton scattering (energy loss) with pT [Q]- We have 
also investigated the effect of isospin asymmetry between p + p and p + p reactions in vr*^ and + 
production and found it to be small. More quantitatively, at = 5.5 TeV the fractional difference 
\daPP - daPP\/daPP varies from 2.5% at pt = 5 GeV to 4.8% at pr = 150 GeV. This is insignificant as 
compared to the projected 50% uncertainty that comes from the extrapolation of lq in LO calcula- 
tions (see Fig. [T]) or the choice of scale in NLO calculations. A recent study showed no deviation from 
DGLAP evolution, Eq. @, at = 10 GeV^ down to x = 10"^ in + iV reactions [Ml] The nuclear 



amplification effect oc A^^^ ~ 10 for a large nucleus is still insufficient to enable measurements of high 
initial gluon density QCD at RHIC, but will play an important role at the LHC. 

2.1 Perturbative QCD hadron composition 

The predicted hadron composition inp + p (p + p) reactions is plotted in the right panel of Fig. ^ The 
proton-i-kaon fraction is seen to increase systematically with px (xt = Ipxl \fs) and is reflected in the 
decreasing 7r'^/0.5(/i+ + /i^). At RHIC and LHC energies this ratio becomes ~ 0.5 at pT — 15 GeV 
and pt — 75 GeV, respectively. At transverse momenta pT — 2 — A GeV the contribution of baryons 
and kaons to + h~ is < 20%. This is significantly smaller compared to data on + reac- 
tions, with the discrepancy being amplified in central A + A. Possible explanations include: enhanced 
baryon production via topological gluon configurations (junctions) and its interplay with jet quench- 



ing ||23|, g4p in ^ + ^ [ ]25| , [2q ], hydrodynamic transverse flow []27[], uncertainty of the fragmentation 
functions Dp/^{zc, Q^) into protons and antiprotons [^], and quark recombination driven by unortho- 



dox (extracted) parton distributions inside nuclei [29]. The approaches in Refs. [M, 26, 27] also address 



the centrality dependence of the baryon/meson ratios in heavy ion collisions at RHIC. In [ pq ] in has been 
shown that similar nuclear enhancement is expected in A, A production (as compared to kaons). The 
combined low-pT baryon enhancement and the growth of the non-pionic hadron fraction in the perturba- 
tive regime may lead to an approximately constant pion to charged hadron ratio in the full measured pT 
region at RHIC at ^/sj^j^ = 200 GeV. We propose that the LHC may play a critical role in resolving the 
mystery of enhanced baryon production in A + A through the significantly larger experimentally acces- 
sible Pt range. Effects associated with baryon transport and transverse flow are not expected to extend 
beyond pj- = 10 — 15 GeV and may result in a detectable minimum of the baryon/meson ratio versus px 
before a secondary subsequent rise. On the other hand, fragmentation functions (possibly enhanced at 
large Zc relative to current parameterizations) are expected to exhibit a much more monotonic behavior. 

3. NUCLEAR MODIFICATION FACTORS 

Dynamical nuclear effects mp + A and A + A reactions are detectable through the nuclear modification 
ratio 

daP^ A daPP in p + A 

dyd^PT dyd^pT 

Rba{pt) = S , , , (5) 

dN^\h) , TAA{b) daPP .^^^^ 

dyd'^PT dyd^pT 

where A and Tj^A{b) = / d^r T4(r)TB(r — b) in terms of nuclear thickness functions Tj^{r) = 
J dzpAiiCjz) are the corresponding Glauber scaling factors [ |30| ] of daPP. We note that in Rba{pt) 
the uncertainty associated with the Kj^io factors, discussed in the previous section, drops out. The ref- 
erence calculations that follow include shadowing/antishadowing/EMC-effect (here referred to as "shad- 
owing"), the Cronin effect, and the non-Abelian energy loss of jets. The scale dependent nuclear PDFs 
read: fa/Ai^iQ'^) = Sa/A{x,Q'^) {Z/A fa/p{x,Q'^) + N/A fa/nix,Q'^)), where we take the isospin 
effects on average and the EKS'98 parameterization ^]\ \ of the shadowing functions S^/Aix, Q^). Initial 



state multiple elastic scatterings have been discussed in ]|3j,|33|,|34j]. From [ ]34j ] the transverse momentum 
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distribution of partons that have undergone an average x = L/\ incoherent interactions in the medium 
can be evaluated exactly for any initial flux dN^^^ (p): 



div(p) = ^ e-^\ / n d\^—^ rfiv(°)(p - qi qn) . (6) 

n=0 '' i=l '^^^ 

Numerical estimates of (^) show that for thin media with a few semi-hard scatterings the induced trans- 
verse momentum broadening exhibits a weak logarithmic enhancement with pT and is proportional to 
L oc A^/^. The transverse momentum transfer per unit length in cold nuclear matter is found to be 
/i^/A ~ 0.05 GeV^/fm [H] from comparison to low energy p + A data [M, 10|. The left top and bottom 



panels of Fig. || show the predicted Cronin-i-shadowing effect m p + Ph collisions at ^/s = 8.8 TeV 
and central Ph + Pb at ^/s = 5.5 TeV without final state medium induced energy loss. The 4% (10%) 
enhancement of Rba at pr — 40 GeV comes from antishadowing and in not related no multiple initial 
state scattering. The observed difference between 7r° and 0.5(/i^ + h~) reflects the different Sa{x, Q"^) 
for quarks and gluons. Cronin effect at the LHC results in slowing down of the decrease of Rba at small 
X as seen in the ^ limit. In contrast, at RHIC one finds ~ 30% enhancement md + Au reactions at 
= 200 GeV and ~ 60% effect in central Au + Au relative to the binary collision scaled p + p result. 
At CERN-SPS energies of y/s = 17 GeV the results are most striking, with values reaching 250% in 
d + Au and 400% in central Au + Au at ~ 4 GeV. For a summary of results on midrapidity Cronin 



effect at the LHC see [g5|l. 

The manifestation of multiple initial state scattering and nuclear shadowing at forward and back- 
ward rapidities y = ±3 in p + P6 at the LHC (for CMS rj < 2.5) and d + Audi RHIC (for BRAHMS 
r/ < 3 ) has also been studied in the framework of a fixed (or slowly varying) initial parton interaction 
strength. At LHC energies at y = +3 (in the direction of the proton beam) the effect of the sequential 
projectile interactions is again small (due to the much flatter rapidity and transverse momentum dis- 
tributions) and is overwhelmed by shadowing, which is found to be a factor of 2-3 times larger than 
the y = result at small px ~ few GeV and vanishes (Rba = 1) at p-r ^ 50 GeV. As previously 
emphasized, initial state gluon showering can significantly change the low-px behavior of the hadronic 
spectra at the LHC beyond the current shadowing parameterization. At RHIC in d + Au reactions at 
= 200 GeV the nuclear modification ratio is qualitatively different. While near nucleus beam (back- 
ward y = —3) rapidity Rba — 0.9 — 1 at forward rapidities y = +3 the nuclear modification factor 
exhibits a much more dramatic px dependence. At small transverse momenta pj^ ~ 1 GeV hadron 
production is suppressed relative to the binary collision scaled p + p result, Rba < 0.8. The maxi- 
mum Cronin enhancement i?^^^ ^ 1.3 (30%) is essentially the same as at midrapidity [Q] but slightly 
shifted to larger pT- We emphasize that both the suppression and enhancement regions are an integral 



part of the Cronin effect [|8), ^, [lOp that is understood in terms of probability conservation and momen- 
tum redistribution resulting from multiple initial state scattering [Q, [16, ^ ^ 34, 35]. At forward 
(in the direction of the deuteron beam) rapidities a calculation as in [^] demonstrates a broader Cronin 
enhancement region with Rba — 25% at p-j- = 5 GeV. This is understood in terms of the significantly 
steeper fall-off of the hadron spectra away from midrapidity that enhances the effect of the otherwise 
similar transverse momentum kicks. While the discussed moderate px interval lies at the very edge of 
BRAHMS acceptance (at y = +3) the same qualitative picture holds aty = +2. 



The full solution for the medium induced gluon radiation off jets produced in a hard collisions 
inside the nuclear medium of length L and computed to all orders in the correlations between the multiple 
scattering centers via the GLV reaction operator approach [36] can be written as (x = w a;/£') 
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Fig. 3: The antishadowing and Cronin effects in p + Pb and central Pb + Pb without energy loss at the LHC 
(y/s = 5.5 and 8.8 TeV) are shown in the left top and bottom panels. The right panel demonstrates the dominance 
of final state radiative energy loss effects at the LHC with a much stronger dependence compared to RHIC. The 
possible restoration of the participant scaling through hydrodynamic-like feedback atpx ^ is also shown 
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where J^l, = is understood. In C 
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are tiie color current propagators, uj 
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tion times, and Az^ = — z^-i are the separations of subsequent scattering centers. The momen- 
tum transfers qj are distributed according to a normalized elastic scattering cross section |wi(qi)p = 
a~i^daei/(Pc[i and the radiative spectrum can be evaluated from for any initial nuclear geometry with 
an arbitrary subsequent dynamical evolution of the matter density. At large jet energies the lowest order 
correlation between the jet production point one of the scatterings that follow has been shown to domi- 
nate and lead to a quadratic mean energy loss dependence on the size of the plasma, AE oc for static 



media [37]. To improve the numerical accuracy for small parton energies we include corrections to third 
order in opacity [Q]. The dynamical expansion of the bulk soft matter is assumed to be of Bjorken type. 



For a summary of results from recent non-Abelian energy loss calculations see |38]. 
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In the Poisson approximation of independent gluon emission [ |39| , |4Q| , |42| ] the probability dis- 
tribution P(e, E) of the fractional energy loss e = J2i ^i/E can be obtained iteratively from the single 
inclusive gluon radiation spectrum dN{x, E)/dx ||40|]. If a fast parton looses e-B of its initial energy 
prior to hadronization its momentum fraction is modified to z* = Ph/Pc{^ ~ ^ ) = ^c/ (1 — e)- The ob- 
servable suppressed hadron differential cross section can be computed from Eq. (Q) with the substitution 



Dh/c{zc,Q^) I deP{e,p,)^Df,/,{z*„Q^) . 

J Zc 



(8) 



The nuclear modification factor Raa{pt) at the LHC is shown on the right panel of Fig. |3| and is com- 
pletely dominated by final state interactions (see left panel). It shows a significantly stronger pT depen- 
dence as compared to RHIC, where jet quenching was predicted to be approximately constant over the 
full measured moderate- to high-transverse momentum range [Q], the result of an interplay of shadowing, 
Cronin effect, and radiative energy loss. The variation of Raa at the LHC is a factor of 5: from 10-20 
fold suppression at pT = 10 GeV to only a factor 2-3 suppression at pT = 100 GeV. The reason for 
such a prominent variation is the hardening of the particle transverse momentum spectra and the insuf- 
ficient balancing action of multiple initial state scatterings. In fact the prediction from Fig. ^ is that the 
suppression in central Ph + Ph at ^/s^vAf = 5.5 TeV at px — 40 GeV is comparable to the factor of 4-5 
suppression currently observed at RHIC. 

The extrapolation of the LHC quenching calculations to small ^ results into suppression 
below participant scaling. More careful examination of the mean energy loss of partons, in particular 
for gluons radiating in nuclear matter at LHC densities, reveals sizable regions of phase space with 
/S.E > E. This indicates complete absorption of jets in nuclear matter. There is experimental evidence 



that this regime of extreme ^naZ state densities may have been achieved at RHIC [ |43| , |44| , |45| ]. In this 
case Eq. (Q) has to be corrected to include the feedback of the radiated gluons into the system. This 
hydrodynamic-like feedback is expected to recover the Npart scaling in the soft px region [Q] - also 
illustrated on the right panel of Fig. |[ The effective initial gluon density derived from the rapidity 
densities used in Fig. | are pg{RHIC) = 30 - SO/fm^ and pg{LHC) = 130 - 275/fm3. These are 
one to two orders of magnitude larger than the density of cold nuclear matter and are suggestive of a 
deconfined QCD state - the quark-gluon plasma. Interestingly, a recent study of non-equilibrium parton 
transport in central Au + Au and Pb + Pb at = 200 GeV and ^/s = 5.5 TeV has found initial 

parton densities corresponding to the lower bound of the intervals quoted above. 



4. CONCLUSIONS 

In summary, a lowest order pQCD analysis of single inclusive hadron production has been performed, 
revealing a systematic decrease with ^/s of the contribution of the next-to-leading corrections to the 
differential cross sections. The predicted da^ / dyd^pr exhibits significant hardening with transverse 
momentum and an increased fractional contribution of kaons and protons at high pT, the latter also being 
true at RHIC energies. In central A + A reactions the nuclear modification factor Raa{pt) at the LHC 
is shown to be completely dominated by final state multi-parton interactions [Q]. For comparison, at 
RHIC Cronin effect and nuclear shadowing also play an important role, leading to an approximately 
constant suppression ratio. At the SPS initial state multiple elastic scatterings dominate, resulting in 
a net enhancement of hadron production. At forward (y = +3) rapidities in d + Au at RHIC the 
Cronin enhancement region is predicted to be broader in comparison to the y = case. In contrast in 
p + Pb at the LHC nuclear shadowing dominates but in order to detect a sizable reduction relative to 
the binary collision scaled p + p cross section measurements at close to proton rapidity (ymax = 9.2 for 
y/s = 8.8 TeV) are needed. 

The predicted decreasing Raa withp^ at the LHC, if confirmed, may have important experimental 
consequences. Comparative large-£^T measurements of the difference in the. full structure of the jet cone 
mp + p and A + A reactions may prove difficult for weak signals and large backgrounds. We emphasize 
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that one of the easiest and most unambiguous approaches for detecting the non-Abelian jet energy loss 
and performing jet- tomographic analysis of the properties of the hot and dense matter created in ultra- 
relativistic heavy ion reactions is through the suppression pattern of leading hadrons. Therefore these 
measurements should enter as an important part of the experimental programs at the LHC. 
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